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Abstract 



J2' 



Gaussian QCD sum-rules are used to analyze all possible two-point correlation functions of scalar gluonic and quark 
currents. The independent predictions of the masses and relative coupling strengths from the different correlators 
are remarkably consistent with a scenario of two scalar states that couple to nearly-maximal mixtures of quark and 
gluonic currents. 
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One of the most fascinating conceptual predictions of 
QCD is the possibility of exotic hadronic states com- 
posed purely of gluons. Such glueball states can have 
the same J PC quantum numbers as conventional quark 
(qq) mesons, and could thus result in supernumerary 
states beyond what is expected from qq structures alone. 
The existence of nineteen scalar J PC = ++ mesons be- 
low 2 GeV (jl is suggestive of an isoscalar glueball in 
addition to a qq nonet and a qqqq nonet. However, the 
observed isoscalars can be mixtures of glueballs and qq 
mesons, which could obscure unique phenomenologi- 
cal signatures of glueballs (e.g., suppressed yy decays). 
It is thus important to determine how the glueball and 
qq content is distributed amongst the observed scalar 
mesons. 

In this paper, we briefly review the Gaussian QCD 
sum-rule analysis of the mixing of scalar gluonium and 
quark mesons ylHl]. Although our findings are comple- 
mentary to QCD Laplace sum-rules [4] (see Ref. ^ for 
a comprehensive review) that conclude that an approxi- 
mately 1 .6 GeV state is a glueball-gg mixture, a unique 
aspect of our approach is a comprehensive combined 
analysis of all possible two-point correlation functions 
that include glueball and qq currentsQ Important in- 
sights that emerge from our analysis include the crucial 
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'Key findings from chiral Lagrangians, lattice QCD, and other 
theoretical approaches provide further support for the scenario of a 
1.5 GeV state with significant gluonium content 1 6, 7]. 



role of chiral violating effects from QCD condensates 
and instantons for a consistent mixing scenario. 

Within a QCD sum-rules approach, the composition 
of a hadronic state may be identified through its cou- 
pling to a field-theoretic current. In particular, we use 
an OZI-inspired definition that qq states \q) couple to 
the light quark current J q = m q qq but decouple from the 
gluonic current J c - aG 2 - aG a G a 



fJV fJV 

(0\J q \q) + , (0\J G \q) * . 



(1) 



Conversely, glueball states \G) couple to the gluonic 
current but decouple from the (light) quark current 



<0|7,|G>*0, <0|/ C |G>*0. 



(2) 



However, mixed states \M) couple to both gluonic and 
quark currents 



<o|/,|m>*o, {o\j g \m)*o. 



(3) 



Three possible correlation functions can then be formed 
from these currents (Q 2 = -q 2 ): 

n cc (e 2 ) = i J d 4 x^-\0\T[J c (x)J c (0)]\0) (4) 
n w (e 2 ) = i j d 4 xe^W[A(x)/ 9 (0)]|0> (5) 

n C9 (fi 2 ) = i j d 4 x^ x (0\T[j G (x)J q (0)]\0). (6) 

In the pure (un-mixed) states scenario there is no in- 
termediate state that couples to both currents, the non- 
diagonal correlator © will be suppressed, and in the 
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absence of an accidental degeneracy, we expect signals 
of different states from the diagonal gluonic correlation 
function (|4]i and the diagonal quark-current correlator 
10. However, if the states are mixtures, then no sup- 
pression of non-diagonal correlator (0 occurs, and sig- 
nals of the same states must occur in all three correla- 
tors. Thus the key signature of mixed glueball-gg states 
is a substantial (i.e., not suppressed) non-diagonal cor- 
relator. 

The original analysis of the non-diagonal correlator 
H has been extended to include contributions of the 
chiral violating mixed condensate and instantons J2l- 
The calculation is quite subtle; one needs to take into 
account the renormalization of the gluonic currenj| 



G\ = [\ + -j ^ j G 2 B - 4^ - {ni u uu + m d ddj B + ... (7) 

and one must also take into account the mixing of con- 
densate coefficients in the operator-product expansion 
|0]. The result for the non-diagonal correlator is H 

n g9 (Q 2 ) =m 2 q Q 2 [A L + A,L 2 ] + m q (qq}C L 

+ m 2 q {aG 2 )^ [B + B\L] + m q (qcrGq) ^ 



Q 2 



- 8 ^ c m q pQ 2 ^Q 2 K X (pQ) K 2 (pQ) 



(8) 



23 (a\ 2 5 (aY (a\ A 
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a a 
D = 4a , B = 6- , B x = -2- , L = log 
n n 



Q 1 



where the dilute instanton liquid model [10] has been 
employed. Eq. ([8]) demonstrates that perturbative and 
gluon condensate contributions have a greater chiral- 
suppression compared with the chiral-violating effects 
of the quark condensate, mixed condensate, and instan- 
ton. Using the definition 



<0|/ G |M> = / C , (0\J q \M)=f q 



(10) 



for the coupling of a mixed state to the currents, the 
(dominant) leading-order perturbative contributions to 
the diagonal correlation functions i4\5\ give Q 



si 



f 



m 2 q E 2 



(11) 



"In this notation, R denotes a renormalized composite operator, B 
denotes bare quantities, and our convention for dimensional regular- 
ization uses D = 4 + 2e spacetime dimensions. 



where £ is a typical hadronic sum-rule scale. Combined 
with ([8]), the approximate contributions to the mixing 
angle 9 from perturbative and gluon-condensate effects 
are chirally-suppressed 



a m„ m„ (aG ) 

sin 26» <s l,sin26> -—r 1 - «: 1 

n E EE 4 



(12) 



However, the quark and mixed-condensate contribu- 
tions to the mixing angle do not experience chiral sup- 
pression 



sin2 ^«<W>, sin2 ^<l^ <<^> 
n E 3 E 5 E 2 £ 3 



(13) 



because the condensate scales (qq) and Mo [see 
Eq. fl% ] are comparable to the hadronic scale E. In 
fact, one can see that the mixed condensate, absent 
in the original analysis HI], is more important than 
the quark condensate. We thus conclude that chiral- 
violating effects can lead to a significant glueball-g^ 
mixing angle, implying that mixing of gluonic and qq 
degrees of freedom has a non-perturbative origin. Qual- 
itatively, this conclusion is similar to that obtained for 
glueball decays |1T] and to that of Ref. lfl2ll which 
demonstrated that instantons can lead to a significant 
mixing between glueballs and (heavy quark) mesons in 
the pseudoscalar channel. 

To go beyond order-of-magnitude estimates of the 
mixing angle, a sum-rule methodology that is sensi- 
tive to multiple states is desirable. Gaussian sum-rules 
(GSRs) are ideal for this purpose. These have the form 
El 



G (s, t) 




exp 



-(t-s) 2 



At 



1 



-p{t) At (14) 



and relate a QCD calculation Go (s, r) to its associated 
hadronic spectral function p(t) with threshold fo. Al- 
though the quantity t (corresponding to the duality in- 
terval) is constrained by QCD, the peak s of the Gaus- 
sian kernel is unconstrained. This implies that as s is 
varied, the resonance peaks of the spectral function are 
probed with equal sensitivity, thereby permitting mul- 
tiple states to be resolved. This behaviour should be 
contrasted with Laplace sum-rules which exponentially- 
suppress excited states, thereby enhancing the ground 
state. 

The original development of GSRs used the heat- 
evolution equation to demonstrate that the finite-energy 



3 Unfortunately, a simple scaling-argument estimate for the instan- 
ton is not possible. 
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sum-rule constraint should be satisfied 11311 . By inte- 
grating both sides of ( TBI with respect to s, it is evident 
that the normalization of Go corresponds to the finite- 
energy sum-rule. Thus the information that is indepen- 
dent of the heat-evolution analysis is encapsulated in the 
normalized Gaussian sum-rules (NGSRs) 



N (S,t) 



Go (s, t) 
M (t) 



M (t) 



/ 



G (S,T)dS, (15) 



related to the spectral function via 



N (s, T, s ) = 



(16) 



A model of two narrow resonances plus QCD con- 
tinuum is used to analyze the NGSR ( fT6b . Although 
superficially this would seem flawed for scalars, reso- 
nance width effects are obscured by the width 2 -\ff of 
the Gaussian kernel for the t ranges where the QCD 
sum-rule is viable. This is illustrated in Fig. Q] which 
compares the contributions to No from a single reso- 
nance of mass m with a square pulse of width T iflill 
centred at the same mass m. It is clear from the figure 
that there is no discernible width effect after integration 
with the Gaussian kernel. This insensitivity to reso- 
nance widths is quite distinct from Laplace sum-rules 
where resonance width effects influence the mass pre- 
diction that emerge from the sum-rules in scalar chan- 
nels |Q. 

N (GeV- 2 ) 




s (GeV ) 

Figure 1 : Comparison of the phenomenological side of the NGSR 
1161 for a single resonance of mass m = 1.0 GeV (solid curves) and a 
square pulse with width T = 0.3 GeV (dashed curves). The overlap- 
ping upper pair of curves is for t = 2 GeV 4 and the bottom (overlap- 
ping) pair corresponds to t = 4 GeV 4 . 

Using the two narrow resonance model, the diagonal 
NGSR (O becomes 



JAtttNo (S, t, so) = n exp - m 2 ) 2 /(4r)] 
+ r 2 exp[-(s - m 2 ) 2 /(4r)] 



(17) 



Since the coupling parameters satisfy r\ + = 1, in 
the diagonal (glueball-glueball and quark-quark) cases, 
where r,- - ff I [f 2 + / 2 2 ), they can be interpreted as the 
relative strength associated with each resonance. 

The resonance parameters and continuum are deter- 
mined by fitting the model to the QCD prediction. The 
QCD expressions for No corresponding to the three cor- 
relators (U-O are too lengthy to be presented here; com- 
plete expressions may be found in Refs. . In addi- 
300 MeV and the PCAC value for the 
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quark condensate, the following QCD parameters are 
used in our analysis (J, 151: 



{aG 2 } = (0.07 + 0.01) GeV 4 , 
(qo-Gq) = M\ (qq) , M 2 = (0.8 ±0.1) GeV 2 , 



(18) 
(19) 



n c = 8.0 x 10~ 4 GeV 4 , p = — GeV 1 , (20) 

0.6 



2.5 MeV < m q (2 GeV) < 5.5 MeV 



(21) 



The results of an independent fit to each correlator 
are given in Table [TJ and shown for the non-diagonal 
(glueball-quark) correlator in Figure [2] The agreement 
between the QCD prediction and the optimized reso- 
nance model is excellent, with no discrepancies that 
would be indicative of additional states or resonance 
width (or resonance shape) effects. The uncertainties 
arising from the QCD input parameters are approxi- 
mately 0.2 GeV for the masses, although there is a cor- 
related effect that leads to a relatively stable splitting 
rri2 - m\ * 0.5 GeV. The uncertainties in the cou- 
pling strengths r, is approximately 0. 1 . Taking into ac- 
count the uncertainties in the mass predictions there is a 
remarkable consistency independently arising from all 
three correlators, providing clear evidence for mixed 
states coupling to both qq and glueball currents, with 
the heavier state coupling more strongly to glueball cur- 
rents. 



Correlator 


m\ 


mi 


r\ 


r% 


so 


JgJc 


0.98 


1.4 


0.28 


0.72 


2.30 


JqJq 


0.97 


1.4 


0.63 


0.37 


2.60 


JcJq 


0.84 


1.4 


0.44 


0.56 


2.75 



Table 1 : Analysis results from the diagonal and non-diagonal NGSRs 
of gluonic and qq currents in the double narrow resonance model. 
Central values of the QCD input parameters have been employed and 
all units are in GeV to the appropriate power. 



The relative strengths r, also reveal a remarkable 
consistency. From the diagonal glueball-glueball and 
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quark-quark results, 



r\ GG) = 0.28 



he r {cc) _ q 72 - he 



f 2 j- f 2 2 
Jig + J2G 

f 2 

jig „(??) 



/lC 2 + /2G 2 



rjw' = 0.63 = — p j , rf ' = 0.37 



f 2 



/l 9 2 + /2 ? 2 



one can determine two possible solutions for the non- 
diagonal glueball-quark case 



JPq) 



fiGflq 



flGflq + flcflq 




(22) 



where the numerical index on the coupling refers to the 
state of mass m, and G, q refers to glueball or quark cur- 
rents. The positive solution is in superb agreement with 
the predicted value for the glueball-quark correlator (see 
final row of Table [TJ. Following Ref. J3l the state cou- 
plings can also be used to determine an effective mixing 
angle <p 

<0|y g |i><0|/,|2> 



tan 



(23) 



<0|/,|2><0|/ 9 |1> 

where |1) and |2) respectively correspond to the states 
with mass m\ and m^. Taking into account the uncer- 
tainties in the Table Q] values leads to <p = 54° ± 4° 
corresponding to states that couple to nearly-maximal 
mixtures of glueball and qq currents. 



N (GeV 




s(GeV-) 



Figure 2: Comparison of the QCD theoretical expression for the 
non-diagonal NGSR N^ q * (S, r, so) with the fitted double narrow reso- 
nance phenomenological model. The upper set of overlapping curves 
are for r = 2 GeV 4 and the bottom set of overlapping curves is for 
T = 4GeV 4 . 

In summary, the three possible correlation functions 
of glueball and qq scalar currents provide signals of 
mixing between scalar glueballs and qq mesons. For 
significant mixing the non-diagonal correlator must be 
non-zero, and this occurs through chiral-violating non- 
perturbative effects from QCD condensates and instan- 
tons. Gaussian sum-rules are able to probe multiple 
states with equal sensitivity, ideal for studying these 
mixing effects. They also show minimal sensitivity to 



resonance width effects, which is important in the scalar 
isoscalar sector where the known states have compara- 
tively large widths. 

The independent analysis of all three normalized 
Gaussian sum-rules reveals a remarkably consistent sce- 
nario of two resonances with a mass splitting of 0.5 GeV 
coupling to nearly-maximal combinations of glueball 
and qq currents, with the heavier state having a slightly 
larger glueball component. This result suggests that 
unique signatures of pure glueball states would be ob- 
scured, complicating definitive experimental identifica- 
tion of glueballs. Our results provide complementary 
supporting evidence for the mixing of qq and gluonium 
to be manifested in the scalar hadronic spectrum as a 
lighter state on the order of 1 GeV and a heavier state 
on the order of 1 .5 GeV 043] . 
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